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I. INTRODUCTION 

Recently, the investigation of semiconductor nanos- 
tructures for photovoltaic applications has been of ever 
growing interest. Potential candidates among these 
low dimensional abs orbe rs are ordered config urat ions of 
quantum wells (QWJP^or quantum dots (QDplil, which 
are widely used in other optoelectronic devices such as 
lasers or light-emitting devices. However, due to the 
unique operating regime of solar cells, optical and trans- 
port properties are equally important, and should not 
be considered in the isolated nanostructure component, 
but for an open system connected to the environment 
via contacts. Since the main attraction of the systems 
is the presence of quantum confinement effects that can 
be exploited to enhance the photovoltaic performance, 
a comprehensive description should be on the level of a 
quantum transport theory. Such a theory was recently 
developed on the example of QW solar cella^*-^ and ap- 
plied to QD solar cells^. 

Excitonic effects play only a minor role in conventional 
inorganic bulk semiconductor solar cells, since in most 
cases, the exciton binding energies are small as com- 
pared to the thermal broadening at room temperature, 
and exciton dissociation is very fast as nothing hinders 
the spatial separation of the carriers. However, this is 
not the case in quantum confined systems, where the 
strong localization of the electron and hole wave func- 
tions leads to a large overlap and thus substantially larger 
exciton binding energies. As a consequence, the ex- 
citonic features in the optoelectronic properties persist 
up to room temperature and have therefore considerable 
impact on the photovoltaic properties of devices based 
on such systems. In the past, excitonic effects in semi- 
conductor nanostructures have been discussed for steady 
state linear absorption or in the regime of high-intensity 
pulse excitation. For the latter, sophisticated quantum- 
kinetic theories were developecPEH. For the description 
of quantum photovoltaic devices, the picture of coher- 
ent excitonic absorption needs to be combined with a 
steady state quantum transport formalism. A suitable 



theoretical framework is provided by the non-equilibrium 
Green's function formalism. However, the shifting of the 
focus from transient interband kinetics to steady-state 
transport does not allow for a straight-forward inclusion 
of excitonic processes: while in the former case coher- 
ent excitonic polarization can be included via the Fock 
term of Coulomb interaction to lowest order, there is no 
equal-time approximation in the latter situation. This 
paper thus aims at the inclusion of excitonic effects into 
a general theory of quantum opto-electronics including 
the transport aspect, which should allow for the study 
of photovoltaic systems where these effects dominate the 
photocurrent response close to the absorption edge. 

The paper is organized as follows. In the section after 
this introduction, the coupling of charge carriers to the 
coherent radiation fields is described based on the NEGF 
theory for a two-band model of a direct gap semiconduc- 
tor, followed by a derivation of the coherent interband 
polarization and the effective interband self-energy due 
to Coulomb-enhanced electron-photon interaction. In a 
further section, these results are used in the simulation 
of the photocurrent respose of bulk and thin film devices, 
where the latter case is represented by a single quantum 
well III-V semiconductor p-i-n diode. 



II. NEGF MODEL OF A CONTACTED 
EXCITONIC ABSORBER 

A. Hamiltonian, Green's functions and self energies 

As a suitable model system, we choose a simple two 
band model of a direct gap semiconductor nanostructure 
selectively connected to ohmic contacts^ and coupled to 
a coherent external photon field, which at this stage is 
treated classically^. Since we are interested in the pho- 
tocurrent response of the system, only the electronic part 
of the system is considered via the Hamiltonian 

H = Ho + H ej + H ep + H ee + He- (1) 
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Hq is the Hamiltonian of the non-interacting isolated 
mesoscopic absorber, H ej is the light-matter coupling, 
H ep encodes the electron-phonon and H ee the carrier- 
carrier interaction. The last term H c describes the (se- 
lective) coupling to contacts required for carrier extrac- 
tion in order to enable photocurrent flow. The charge 
carriers in the two bands b = c,v are described by the 
field operators ^(r, t), which provides the Hamiltonian 
representation 

«(<)= E / d 3 r& a (r,t)H(r,t)i> b (r,t). (2) 



a,o— c,v 



The renormalizing effects of the interaction and con- 
tact Hamiltonian terms on the isolated system are ex- 
presse d wit hin non-equilibrium many-body perturbation 
theory^ 11221 in terms of corresponding self-energies en- 
tering the generalized Kadanoff-Baym equations for the 
charge carrier non-equilibrium Green's functions (1 = 
{ri,ii e C},C: Keldysh contour in the complex plane, 

f di = J dPnfcdh) 

Gq 1)G(1, ]/) =5(1,11) 

+ J d2 £(1,2)0(2,10, (3) 

[Gj]" 1 ^, V)G(h20 =8(1,20 

+ J d2G(l,2)E(2,lO, (4) 

where 

[GoHLlDU = Nfi^ - [ffo(ri)U S(l,r_)S ab (5) 
and the contour-ordered Green's functions are defined via 



Gab{ht) 



(*a(l)^(l'))c 



(G) 



for band indices a, b. The self-energy term in the above 
equations for the Green's functions may be divided into 
the contributions from the interactions and the contact 
term, 



£(l,2) = £ 7 (l,2) + 5] c (l,2), 



(7) 



where the interaction term contains the effects of 
electron-photon, electron-phonon and electron-electron 
coupling, 

S J (1,2) = S e7 (l, 2) + E ep (l,2) + E ee (l,2). (8) 

Following the standard real-time decomposition rules^ 
applied to ([3| and a special band decoupling procedure 
described in App. [Aj the equations for the retarded, 
advanced, lesser and greater components of the non- 
equilibrium Green's functions for charge carriers can be 
written in the standard intraband form used in transport 



calculations, (1 = {i"i,ti e R}, / dl = J d?r 1 dt\) 
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Go* (1,2)-£*(1,2) G« (2, 10 = 5(1,1') 



G<(1,1')= J &2j d3G*(l,2)E<(2,3)G^(3,l0, 

a e{c,v} (9) 



where 



^(1,10 =E i 00 (l,l') + E£(M')> 

i = B,, A, «s, a G {c, v}, 



(10) 



with the effective band-coupling self-energy from the sin- 
gular contributions given by 

££(l,l')=y d2 J d3Ei b (l,2)G 4 bb (2,3)Ef a (3,l'), 

i = R,A,^, a ^ 6, a,be{c,v}. (11) 

In the situation under investigation, the singular inter- 
band self-energy itself is of the form 

£ 5 a6 (l,2) = [E^n.txWn -r 2 ) + £*f (r ls r a ,tO' 

x<5(ti-t 2 ), (12) 

where E^' fc e 7 encodes the coupling of electrons to a coher- 
ent photon field and E^' & ee is the (non-retarded) Fock term 
of the (screened) Hartree-Fock approximation to carrier- 
carrier interaction leading to Coulomb enhancement of 
the optical interband transitions due to the electron-hole 
coupling. 



B. Singular self-energy and coherent interband 
polarization 

The self-energy due to the light-matter interaction can 
be written in terms of the vector potential of the classical 
electromagnetic field, 



(r,t) 



to 



(A(r,t))-p(r) 



(13) 



with (A(r, t)) — A co /j(r,t). The Coulomb term is 



E*nn, r 2 , t) = ihV{v x - r 2 )G< (n, r 2 , t, t + ), (14) 

where V is the (screened) Coulomb potential, and 
depends thus on the coherent interband polarization 
through the interband Green's function, which the de- 
coupling provides in the form 



G<(1,1') 



d2 / c23 



G«(1,2)E* C (2,3)G<(3,1') 



G<(1,2)E^ C (2,3)G^(3,1') 



- / d2 I d3E* c (2,3)r M (l,2,3,l') 



(15) 
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where we have defined 

r„ c (l, 2, 3, 1') =ift[<7* (1, 2)G< (3, l') 
+ Gw(l, 2)G^,(3, 1') 



(17) 



Under steady state conditions, Fourier-transform to the 
energy domain yields 



E* (r 1 ,r a ;B)= / d 3 r 2 / d 3 r 3 



2tt/i 

i = ii, A, % (18) 



with the singular self-energies given by corresponding 
Fourier transforms of Eqs. (13 1 and (14). Inserting 



the explicit expressions for the latter leads to a Bcthc- 
Salpeter type equation (BSE) for the coherent polariza- 
tion, which for steady state reads 



G< c (t u t v ;E) 



d 3 r 2 / d 3 r 3 



A coh (r 2 , E) • p(r 2 )<S(r 2 - r 3 ) + ihV(r 2 - r 3 )G< (r 2 , r 3 , E) 

x % c {r i,r 2 , r 3 , t v ;E) (19) 
= G< c{0) (r 1 ,v 1 ,;E) + J d 3 r 2 J d 3 r 3 V(r 2 - r 3 )% c (r u r 2 , r 3 , r v ; E)G< (r 2 , r 3 , E), (20) 



with 



d 3 r 2 A coh (r 2 ,E) ■ p(r 2 ) 



m h 

xT vc (r 1 ,r 2 ,r 2 ,r 1 ,;E) i21) 

i 

~~~h, 

x 5(r 2 - r 3 )% c (T 1 ,T2,T3,r v ;E) 

(22) 



J d 3 r 2 J d 3 r 3 ^(y 2 ;E) 



the coherent interband polarization of non-interacting 
electron-hole pairs, where 



7; c (ri,r 2 ,r 3 ,r 4 ;£;) = ih 



dE' 
2^h 



G*(T 1 ,r 2 ;E , -E) 



x G< (r 3) r 4 ; E') + G< (n, r 2 ; £' - £)G^(r 3 , r 4 ; £') 

(23) 

Here, it is interesting to note that T vc {*i, r 2, r 2 , r 4 ; £J) = 
T^J (ri , r 2 ; E) is the retarded component of the random- 
phase approximation of the incoherent polarization func- 
tion used to describe the interband coupling that is not 
singular in time^. The microscopic, non-local interband 
susceptibility \vc is introduced via 



G< c (r 1 ,ry,E)=- / d 3 r 2 / d 3 r 3 d vc (r 2 - r 3 ) • £(r 2 , E) 



x Xvc(ri,r 2 ,r 3 ,r v ;E), 



(24) 



where d = — er is the dipole operator. The linear 
macroscopic interband susceptibility is obtained from the 



r 



macroscopic interband polarization given b 



P(r l ,E)=-ih I d 3 r v \d cv {v 1 ^v v )G< c {v 1 ,Y V -E) 

+ d„ c (r 1 -r r )G<(r 1 ,r r ;£)] (25) 
= J d 3 ry^(ryr y ;E)£(ry,E), (26) 



i.e. 



<b-s- 



5?(n,r 1 /;B)= / d 3 r 2 J d 3 r 3 |d CT (n - r 2 )d vc (r v - r 3 ) 
x Xvc(ri,r v ,r 3 ,r 2 ;E) + d vc (r 1 - r 2 )d cv (r v - r 3 ) 



x x ct ,(ri,iv,r 3 ,r 2 ;£0 



(27) 



The susceptibility can be used to compute the local ab- 
sorption coefficient, 



E, 



1 o 



az(r;E 7 ) = _ 



£oc 0y /ebE 1 \m Q J r 



$Xii(r,r;.E 7 ), i € {x,y,z}, (28) 

2 



lim p,-(r)&(i/)G;f c (r,r';.EL), 



(29) 



where Co is the speed of light in vacuum and e& is the 
background dielectric constant. The corresponding av- 
erage (bulk) absorption coefficient may then be defined 
via «i(£ ; 7 ) = V -1 J d 3 raj(r; -E 7 ^, with V the absorbing 
volume. 
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A BSE-type equation similar to (20) can also be de- with 



rived for the singular self-energy using (16) in Eq. (14) 
E^Cn , r'j , E) =£* 6>(0) ( ri , ri ,E) + V(r 1 - r[) 



x y" rf3r 2 J d 3 r 3 Tab(ri,r 2: r 3 ,r 1 r 1 E) 
xZ s ab (r 2 ,r 3 ,E), a^be{c,v}, (30) 



where S fb,(0) 



With the carrier Green's 



function modified by the effective band coupling self- 
energy via Eq. (9j), the spectral response S(hu) 7 ) = 
J (fk/J-y) /[e0 7 (/iu; 7 )]7 where <^> 7 denotes the spectral pho- 
ton flux, is determined from the steady state current in- 
duced under monochromatic illumination in the interact- 
ing region, for which the standard Meir-Wingreen expres- 
sion is usecP3 



s,, 



dE 
2^ 



tr^iry-E) 



x {n F (E - /i a )A(r', r; E) + iG< (r', r; E)} 



(31) 



where S a is the surface area, T a the broadening function 
and Ha the chemical potential of contact a, np is the 
Fermi function and A = i{G R — G A ) the spectral function 
of the fully interacting and contacted absorber, which 
here may be either bulk or a thin film of semiconductor 
material. 



III. APPLICATIONS 

In this section, the self-consistent band coupling self- 
energy approach derived above is first validated for the 
case of a contacted bulk absorber and then implemented 
in the existing NEGF model for thin-film and quantum 
well solar cells^H 



A. Bulk 

For a periodic bulk material, the BSE for the interband 
polarization function can be rewritten in Fourier space as 
follows: 

G<(k;£) = G< i(0) (k;£) 

+ T(k;£)^V(k-q)G<(q;£) (32) 



G< cm (k;E)-^Tik;E)^ b (k,Ei (:>,:!) 



: T(k;E)X s vc (k,E). 



(34) 



If the photon momentum is neglected as compared to the 
electron quasi-momentum, the non-interacting polariza- 
tion function reads 

G< j(0) (k; E) -^Yl ^coh(<L E) ■ p„ c T(k; E) (35) 



T(k;E) =ih 



dE' 



G*(k;E'-E)G<(k;E') 



2nh 

G< v (k;E'-E)G A s (k;E') 



(36) 



To lowest order, inserting the quasi-equilibrium approx- 
imation for the bulk Green's functions, 



G< a (k;E) =2mf(E)S(E-s a (k)), 
G*l A (k-E)=[E-e a {k)± lv ]-\ , 
the following expression is obtained 

/„(k)-/ c (k) 



T(k;E) = 



e c {k)-E v {k)-E + ir}' 



(37) 
(38) 

(39) 



which, used in (35), leads to the standard form of the 
macroscopic polarization function^. 

Assuming complete isotropy, one may neglect the an- 
gular dependence and arrive at the equation 

G<(k;E) = G< cm (k;E)+T(k; E) J dqV(k,q)G< c (q;E) 

(40) 



with the effective Coulomb potential 



V(k,q) 



47T 2 ££ k 



k 2 



2kq 



k 2 + q 2 - 2kq 
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(41) 



The BSE in (40) can then be rewritten as 

G&, (0 )(fc;£) = / dqM(k,q)G<(q;E), (42) 

with 



M(k,q) = S(k - q) - T(k;E)V(k,q), 



(43) 



which can be solved in discrete momentum space via in- 
version of matrix A4. Fig.JT^,) shows the bulk absorption 
coefficient 



a{E-y) 



e \ -2 

m J 



^G<(k ;J B 7 ) (44) 



for a two band effective mass model of a direct semicon- 
ductor with and without Coulomb correlations as derived 
via the macroscopic susceptibility from the coherent po- 
larization function given in ( |26[ ). The parameters used in 
the simulation are given in Tab. [l]and with exception of 
the hole mass correspond to GaAs. The inverse screening 
length is taken at qo = 10 6 to _1 . 

The expression for the (singular) interband self-energy 
may now be rewritten using the above result ( 33 1 for the 
interband Green's function, 

£* c (k, E) =SS'r(k, E) + Ejf b (k, E) 



— y° 



,(0) 



.(0) 



(k,E) + ihJ2VQi-c l )G<(< l ;E) 
q 

(k, E) + V(k - q)T(q; S)S* c (q, E), 

(45) 
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In order to obtain the Coulomb enhancement factor 
for the effective interband coupling, the equation is 
formulated for the normalized self-energy <r* c (k, E) = 
[£* c , ,(25)] -1 £* c (k, E), neglecting the quasi- momentum 
dependence of the momentum matrix elements, 

a 5 vc (k, E) =1 + V(k - q)T(q; S)^ c (q, E), (46) 
q 

which is idependent of the exciting field an hence related 
to the macroscopic interband susceptibility. The effec- 
tive interband self-energy for monochromatic illumina- 
tion with frequency loq can then be written as follows: 

£<^(k; E) H£^ (0) (^ )| Vi(k, ^o)| 2 

xG<(k;E-huj ). (47) 

The spectral response displayed in Fig. [TJd) is obtained 
from the photocurrent given by the bulk version of (31 1, 



J 7 (fiw 7 ) = 



x A(k;E) + iG<(k;E)} 



(48) 



where the Coulomb corrections enter via the effective in- 



terband self-energy ( 18 ) used in the equations for the 



Green's functions G°, a 



A2 
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- 
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FIG. 1. (a) Linear absorption coefficient of a bulk direct gap 
semiconductor absorber described by two parabolic bands, 
as computed via the coherent interband polarization func- 
tion corrected for electron-hole interactions. For comparison, 
the absorption of non-interacting electron-hole pairs is shown 
as well, (b) Spectral response for the same system as com- 
puted via the effective electron-photon self energy corrected 
for electron-hole interactions. 



TABLE I. Material parameters used in the simulations 





m*/mo 


ml/mo 


E g 


Pcv/mo 


e 


GaAs 
ALGai-^-As 


0.067 
0.095 


0.1 
0.1 


1.42 eV 
1.82 eV 


18 eV 
18 eV 


13.6 
12.2 



B. Thin films 

If periodicity is restricted to the transverse dimensions, 
Eq. (20) becomes 



G< c (k\\, zi, z[; E) = G< c (0) (kj|,zi,zi;£:) 
dz 2 J dz 3 T{\\,z ll z 2l z 3 ,zy 1 E) 

x Y J V ( k \\ -CL h z 2 ,z 3 )G< c (q h z 2 ,z 3 ;E), (49) 
ill 

where V is again the screened Coulomb potential. With 
the approximation of angular isotropy in the transverse 
dimensions, the BSE equation may be written 

Ctfc(^|| > z i; z i> E) — ^dc,(o) (fyl i z ii z 'iiE) 
dz 2 J dz 3 T{k\\ , zi,z 2 ,z 3 , z v ; E) 

x J dq\\V{k\\,q\\,z 2l z 3 )G< c {q\\,z 2 ,z 3 ]E), (50) 

with the effective, statically screened Coulomb potential 
=2 r^TT e -g(fc||,<?n,0)|*i-* 2 | 



87T 2 ££ 



V(k h q h z 1 ,z 2 ) 

q(k .,q .0) =Jk 2 



1\\ 



d6 



2h 



cost 



9o 2 - 



(51) 



Neglecting the short-range contribution of the Bloch 
functions, the Coulomb-matrix element for a localized 
real-space basis set may be approximated as Fy (fen , gy) « 
F(fc||, qu , Zi, Zj). Fig. 2 shows the spatial and momentum 
dependence of the effective, statically screened Coulomb 
potential for two different values of the inverse screening 
length, go = 10 6 m^ 1 and go = 10 9 m _1 , corresponding 
to the limiting cases of weak and strong screening, re- 
spectively. As is to be expected, a small screening length 
leads to a long-range interaction that is strongly localized 
in momentum space, with the opposite behavior in the 
case of large screening length. 

Using the localized basis representation, the BSE for 
the coherent interband polarization function may be ex- 
pressed as a matrix equation in analogy to the bulk case, 



[G< c (E)] a = [G< (E)] C 



^2M aP (E)[G<(E)]p, (52) 



where multi-index notation is used, with a = fen) 
and (3 = (I, m, qu), and 

M a{ijk[[ ),f3(i mqil )(E) =Tii m j(kf,E)V lm (k\\,qi\)Aq\\. (53) 

The corresponding equation for the linear susceptibility 
X may be obtained from the above equation 

[Xvc{E)] a = [ Xvc , {0) iE)} a +J2M a p(E)[Xvc{E)}p, (54) 





(> 



with XvcJO) = T ■ This leads to the equation 

~l-M(E)] Xvc (E)=x vcA0 )(E), (55) 



where 



Rot (-E0] a (y fey) = y]Xvc,aij{k\\', E). 



(56) 



The expression corresponding to ( 45 1 for the self- 



consistent equation for the singular self-energy reads 
E* c (k||,si,si;.E) = E* C)(Q) (ky,zi,4;£7) 

+ 5^V(k||-q||,2!i-a4) 



x J dz2 J dz 3T(q\\,z 1 ,z 2 ,z 3 ,z[;E) 

xE s vc (d lb z 2 ,z 3 ;E). (57) 

In the discrete basis, the angular isotropy limit of the 
above equation is 

Imqu 

xTu mj (q\\' E )]^ S vc}lm(q\\;E), (58) 

which can again be solved via matrix inversion. 

In general, the size of the matrices appearing in Eqs. 



( 52 1 and ( 58 1 prohibits the computation of the full ma- 



trix. In a first approximation, off-diagonal elements in 
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FIG. 2. (Color online) Quasi-momentum [z = z' , left col- 
umn) and spatial (qu — k\\ , right column) dependence of the 
statically screened Coulomb interaction for real space resolu- 
tion in one dimension. In the low-screening limit (qo = 10 6 
m -1 , upper row), the interaction is local in momentum space 
and uniform in real space. In the strong screening regime 
(go = 10 9 m _1 , lower row), the situation is reversed. 
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FIG. 3. (Color online) a) Coherent excitonic absorption and 
spectral response of a 5 nm GaAs quantum well embedded 
in the intrinsic region of a Al^Gai-^As pin diode, b) Local 
density of states of the same system, revealing the situation 
at finite field and substantial scattering. 



the spatial indices (i, j) are neglected, the long-range con- 
tributions of the electron-hole Coulomb interaction are 
thus lost, which may lead to an underestimation of the 
exciton binding energy. The approximation is reasonable 
in the case of strong screening, especially for the polar- 
ization function, where a local interaction removes the 
second spatial integration. In the equation for the singu- 
lar self-energy, only the diagonal elements are modified 
by the Coulomb interaction if the potential is local. In 
the other extreme of low screening, where the Coulomb 
potential is spatially constant, one still needs to account 
for the non-locality of T ■ In the following numerical ex- 
amples, the spatial integration over z 3 in (49 1 is consid- 



ered via replacing the interaction potential with the term 
V(k\\,q\\;z 2 ) = / dz 3 V(k\\ , q\i ; z 2 , z 3 ). For consistency, 
the same correction factor is used for the effective po- 
tential in ( 57 ) . Fig. |3^) shows the effects of the Coulomb 
correlations on absorption and photocurrent response of 
a 5 nm wide GaAs quantum well embedded in the center 
of the intrinsic region of an Al^Gai^^As (x ~ 0.3) p-i-n 
diode at a contact Fermi level splitting of 1 V and for 
go = 10 6 m -1 . The parameters for the bulk materials 
are given in Tab. [TJ and the band offsets are AE C = 0.2 
eV and AE V = 0.15. The band bending is obtained from 
self-consistent coupling to Poisson's equation. Scattering 
is treated as in Ref.El The LDOS of the quantum well re- 
gion displayed in Fig.^) reflects the effects of the sizable 
built-in field of ~ 100 kV/cm and the strong electron- 
phonon interaction. The corrected susceptibility shows a 
distinct exciton peak broadened by phonons. The spa- 
tially resolved susceptibility in the QW, proportional to 
the local generation rate, is shown for free-electron hole 
pairs in Fig. [4^,) and for excitons in Fig. [4Jd). Again, 
both the appearance of the exciton peaks below the ab- 
sorption edges of the non-interacting system as well as 
the enhancement of the quasi-continuum absorption rep- 
resent the salient features of the correlated transitions. 
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FIG. 4. (Color online) a) Coherent local susceptibility in the 
QW region for the case of non-interacting electron-hole pairs, 
b) Same for the system with excitonic corrections, exhibit- 
ing the salient features of correlated transtions, with exciton 
peaks broadenend by electron-phonon interaction and an en- 
hancement of the quasi-continuum. 



IV. CONCLUSIONS 

In this paper, a consistent inclusion of excitonic ef- 
fects into the computation of the photocurrent response 
of photovoltaic nanostructures was presented. While a 
full treatment of the two particle interactions is still out 
of range in the context of quantum transport simulations, 
the excitonic enhancement of the coupling to classical 
radiation fields can be considered via the correspond- 
ing modification of the electron-photon self-energy enter- 
ing the equations for the charge carrier non-equilibrium 
Green's functions. However, since the correlations near 
the band gap do also have a strong impact on any in- 
terband recombination process, it remains desirable to 
extend the theoretical treatment to incoherent excitons 
formed by electronically or optically injected carriers. 
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Appendix A: Band decoupling procedure 

The standard real-time decomposition rules^ applied 
to ([3| yield the coupled equations for the retarded com- 
ponents of the intra- and interband Green's functions, 

Go,L(l, ^(l, 1') =8(1, 1') + / d2 [E« (1, 2)G* (2, 1') 



E«(1,2)G«(2,1') 



GnL(M)G«(l,l') 



(12 



E* c (l,2)G«(2,l') 



+ £«(1,2)G&(2,1') 
Introducing the new quantity 

fiR _ _ yR 1 _1 



(Al) 



(A2) 



(A3) 



in ( A2 ) , the retarded interband Green's function can be 



written as 

iR f-\ t'\ — I An I ai nR fi o*\v<5 



G*(l,f) = / d2 / d3G«(l,2)£* c (2,3)G£(3,l'). 



(A4) 



Inserting the above expression in (All provides a closed 



equation for the retarded intraband Green's function, 

Gl (1, 1') = [G -L(l, 1') - Ef c (l, 1') - Eff (1, 1')] 1 

(A5) 

= G -i c (l,l')-Ef c (l,l')]~\ (A6) 

where the contribution from the singular terms to effec- 
tive band-coupling intra-band self-energy E is 

E^(l,l') = J 62 J d3E^(l,2)G«(2,3)E* c (3,l / ). 

(A7) 

In the same way, the lesser and greater components of 
the Green's functions can be decoupled: starting from 



Go,L(li 1)G< (1, 1') = J d2 [Ef^l, 2)G< (2, 1') + Ef c (l, 2)G< (2, 1') + E< (1, 2)G£(2, 1')] , (A8) 
G :L(M)G<(1,1') = / d2 [Et(l,2)G<(2,l') + E«(l,2)G<(2,l') + E<(l,2)G4(2,lO] . (A9) 



the interband correlation or coherent polarization func- 
tion is written as 



G^ c (l, 1) 



d2 / d3 



<?£(!, 2)E* C (2,3)G< (3,1') 



where 
G<(1,1') 



d2 / d3G«(l,2)E<(2,3)G^(3,l') 



(All) 



■G<(1,2)££ C (2,3)G£(3,1') 



(A10) 



was introduced. Replacing the interband term in (A8) 
then yields the intraband correlation function 



with 



G<(1, 1') 



d2 J d3G£(l,2)[E<(2,3) 

Ef<(2,3)]G c A c (3,l') (A12) 

d2 J d3G*(l,2)E< (2, 3)0^(3,1') 

(A13) 



^ 5 <(1,1')= / d2 / d3 £*,(!, 2)G< (2,3^(3,1'). 



(A14) 



The expressions for the valence band self-energy correc- 
tions are obtained from analogous derivations as and are 
identical to the above result with c -H- v. 
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